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The function E — F(v) expresses the dependence of a discrete eigenvalue E of the Schrodinger 
Hamiltonian H = —A + vf(r) on the coupling parameter v. We use envelope theory to generate a 
functional sequence {/' fe '(V)} to reconstruct f(r) from F(v) starting from a seed potential /'°'(V). 
In the power-law or log cases the inversion can be effected analytically and is complete in just 
two steps. In other cases convergence is observed numerically. To provide concrete illustrations 
of the inversion method it is first applied to the Hulthen potential, and it is then used to invert 
spectral data generated by singular potentials with shapes of the form f(r) — —a/r + bsgn(q)r q 
and f(r) — —a/r + &ln(r), a, b > 0. For the class of attractive central potentials with shapes 
■ f(r) = g{r)/r, with g(0) < and g'(r) > 0, we prove that the ground-state energy curve F(v) 

determines f(r) uniquely. 
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I. INTRODUCTION 

The purpose of this paper is to report a robust method for reconstructing the detailed shape of a singular attractive 
potential from some of the bound-state spectral data that it generates. The general setting for the problem is the 
^ I discrete spectrum of a Schrodinger Hamiltonian operator 
ON 

H = -A + vf(r), r=||r||, (1) 

where f(r) is the shape of an attractive central potential, and v > is a coupling parameter. We shall assume that the 
potential is monotone non-decreasing and no more singular than the Coulomb potential f(r) = — 1/r. The arguments 
we use apply generally to the problem in d > 1 spatial dimensions, but, for dehniteness, we shall usually assume that 
d = 3. The operator inequality [l], [1[ 

-A> d>3, (2) 



- 1—1 

X 

implies that a discrete spectrum exists for sufficiently large coupling v > 0. For d — 3, the Hamiltonian H is bounded 



below by 



E > min 

r>0 



(3) 



and a simple trial function can be used to establish an upper bound. Thus we may assume, in particular, that the 
ground-state energy may be written as a function E = F(v). An explicit example of the class of problems we consider 
is provided by the Hulthen potential whose shape is given by /(r) = — l/(e r — 1) and whose s-state (£ = 0) eigenvalues 
E n are given Q exactly for d = 3 by the formula 

( v — v? \ 

E n =F n (v) = - ( , v>n 2 , n = 1,2,3,.... (4) 
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The problem discussed in the present paper may be stated as follows: given, for example, the curve F\(v), can we 
use this spectral data to reconstruct the potential shape f(r)7 We call this 'geometric spectral inversion'. It must 
be distinguished from 'inverse scattering theory' (34l^|. and, in particular, from 'the inverse problem in the coupling 
parameter' Q. In this latter problem, the energy eigenvalue E is fixed, and the data to be inverted comprises the 
set of all possible couplings leading to a radial eigenstate with this given energy, along with certain scattering data. 
In geometric spectral inversion, we simply use one spectral curve, say F\(v). Progress has been reported on this 
problem in a sequence of earlier papers [13l4l7j |. From the point of view of problems in d = 1 dimension, we were 
at first concerned that for excited states, which often only exist for v sufficiently large, some essential pre-image 
data might be 'missing'. However, the analytical inversion of the WKB approximation [15| showed that the WKB 
inversion process improves in quality as one goes to higher states, and is moreover asymptotically exact. Progress has 
been made on a related problem, namely on inverting spectral data of the form E(£), where £ is the orbital angular 



momentum quantum number, for fixed coupling v\ this work by Grosse and Martin is described in Ref. 18|. 

We suggest two areas of application for the form of spectral inversion F{v) — > f(r) which is the main topic of this 
paper. The spectrum generated by the outermost electron in a large atom is often estimated by use of a screened- 
Coulomb potential with parameters dependent on the atomic number Z . It is envisaged that atomic spectral data 
for a sequence of atomic numbers could be used to reconstruct an appropriate effective potential for the model. In 
another area, quantum field theory can be used to predict the dependence E — F(v) of the energy of a two-body 
bound system on the coupling v. Geometric spectral inversion could then be used to construct from F(v) a potential 
f(r) that would reproduce the same spectrum in a quantum- mechanical model. The idea is that instead of fitting the 
parameters in an appropriate potential model, we instead develop an inversion method that allows us to reconstruct 
the potential shape directly from an eigenvalue curve. 

In the present paper we study the inversion of spectral data arising from singular potentials. We develop a functional 
inversion algorithm based on envelope theory, a geometrical theory of energy functions such as F n (v). This was first 
introduced }19| in 1980 as a method for finding energy bounds for many-body problems, and has since led to a 
number of refinements and applications (20H2^ | , including semirelativistic problems [25M29j | . In sections II and III we 
first review enough of this theory to serve our present purpose and to make this paper essentially self-contained. In 
section IV we develop the inversion algorithm. We then test the algorithm on some concrete problems. In section V 
we show that the algorithm inverts spectral data arising from pure powers, or the log potential, exactly, in just two 
steps. In section VI we invert the Hulthen spectral data from Eq. Q. In section VII we apply the inversion algorithm 
to the spectral data generated by the family of potentials f(r) = —a/r + bw(r), where w(r) — r, w(r) = r 2 , and 
w(r) — ln(r). In these latter problems, for which analytic expressions for the spectral curves F(v) are not available, 
we begin in each case with the potential /(r), find F(v) numerically (by shooting methods), and then apply the 
sequential inversion method presented in this paper to reconstruct f(r) again from F(v). The spectral inversion 
F — > f is the principal concern of the present work. In section VIII we consider the ground state energy curve F(v) 
for the problem H = —A + vg(r)/r, where g(0) < 0, and g'(r) > 0. For this class of problems we prove that F(v) 
determines f(r) = g(r)/r uniquely. 



II. EXACT REPRESENTATION OF SPECTRAL FUNCTIONS BY KINETIC POTENTIALS 

The discrete spectra of operators such as H = —A + vf(r), which are bounded below, may be characterized 
variationally [3C|. Thus, the ground-state energy may be written 

F(v) = inf (ib, Hip). (5) 

V>£T>(ff) 

IM>ll=i 

Since H depends on the coupling v, so therefore does the domain V(H). However, for the problems considered, either 
H has discrete eigenvalues, perhaps for v greater than some critical coupling vx, or the entire spectrum of H is discrete 
for v > 0. The kinetic potential f(s) associated with the potential shape f(r) is defined (for the ground state) by a 
constrained minimization in which the mean kinetic energy s = (—A) is kept constant: 

f(s)= inf (1PJ1P). (6) 

(i)>,— Ai/))=s 

The eigenvalue F(v) of H is then recovered from f(s) by a final minimization over s: 

F(v) = mm [s + vf(s)] . (7) 
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The spectral function F(v) is concave (F"(v) < 0); moreover, we have shown [l3[ that 



F"(v)f"(s) = 



(8) 



Thus F(v) and /(s) have opposite convexities and are related by the following Legendre transformations / -H- F [31| 



f(s) = F'(v), s = F(v)-vF'(v), 
l/v = -f'(s), F(v)/v = f(s) - sf'(s). 



(9) 
(10) 



F(v) is not necessarily monotone, but /(s) is monotone decreasing. Eq. © enables us also to work with the coupling 
as a minimization parameter. For this purpose we write the coupling as u and we have from Eq. ([7]) 

F(v) = min \F(u) - uF'(u) + vF'fu)] . (11) 

This is particularly useful in cases where f(s) is difficult to find explicitly. 

By considering finite-dimensional linear spaces in T>{H) we can extend these definitions and transformations [HIHH 
to apply to the excited states. For example, in d = 3 dimensions, if D n is an n-dimensional linear space of radial 
functions {&} contained in T>(H) and in the angular-momentum space labelled by Y™, then we have 



Ent = inf sup (ip,Hip). 

II V> 11=1 



(12) 



We now scale the linear space D n so that we can fix (—A) = s. We define a scaling operator a by (aip)(r) = tp(ar), a > 
Then, if we let 



we may then define 



&D n = span{(T^j}^ =1 , 
9 n = |J oD n . 



(13) 



(14) 



er>0 



We note that this union of linear spaces is not itself a linear space. We may now define the excited-state kinetic 
potentials by 



fru(s)=w£ sup (ipjip)- 
IHI=i 

(•0,-A«)=s 



(15) 



We shall not usually need to use this abstract definition since Eq. © enables us to generate kinetic potentials directly 
from known energy functions F(v). For example, with the s-states for the Hulthen potential /(r) = — l/(e r — 1), we 
have from Eq. (|4|) 



fn(s) 





'As 


i 

2 




— + 1 


- 1 




n z 



n = 1,2,3,.. 



(16) 



The point of these alternative expressions for the spectral curves will become clear when we discuss approximations 
in the next section. Another form of expression, useful for our present task, is obtained if we change the kinetic-energy 
parameter from s to r itself by inverting the (monotone) function f(s) to define the associated if-function by 



KW(r) = s=(f- 1 °f) (r). 



(17) 



Now the energy formula Eq. (J7J becomes 



F(v) = min \KU\t) +vf(r) 

r>0 L 



(18) 



A sleight of hand may be perceived here since K depends on /. However, we do now have a relation that has F on 
one side and / on the other: our goal is to invert this expression, to effect F — > /. We shall do this in section (4) by 
using a sequence of approximate K- functions which do not depend on /. 
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Another class of soluble problems is the set of pure-power potentials the form of whose energy functions are 
determined by scaling arguments. We have [2l| 



f(r) =sgn(< Z )r« = E rU (q)v^, q > -2, q + 0, 

where the unit-coupling eigenvalues E n i(q) are known in special cases, such as the Coulomb problem E n i 
-l/(4(n + if) and the harmonic oscillator E nl (2) = in + 21 - 1, where n = 1, 2, 3, . . . and ^ = 0, 1, 2, . 
corresponding kinetic potentials for this family become from Eqs. ([9]) and (fl9|) 



(19) 

-1) - 
The 



qE ne (q) 



f(s) = - 

q 



By contrast, the K- functions are much simpler and are given by 

K^(r) = 

where 

2 



Pnt(q) = ^ 



For example 



p nt (-l)=n + £ and P„ £ (2) = 2n + ^ 



(20) 



(21) 



(22) 



(23) 



By defining P n i(0) = lim^o Pni(q), we not only make P n i(q) continuous at q = but we also exactly accommodate 
the ln(r) potential [2l|. It has been proved that the P(q) functions are monotone increasing |22| . and they appear to 
be concave. By contrast, the functions E(q) are more complicated and have infinite slopes at q — 0. If we denote the 
eigenvalues of —A + ln(r) by E^, we have by scaling 



/(r) = ln(r) F ni (v) = --vhx(v/v nt ), v nl = exp(2£^). 
Meanwhile, the corresponding kinetic potential and i-T-function are given by 

/»=2&-i]n(2e a ) and K^(r) = 



(24) 



(25) 



Numerical values for P n g (q) may be found for d = 3 in Ref. [32[ and for more general d > 1 in Ref. [33[ . 

The remaining task for this section is to note the scaling and invariance properties of f(s) and K^\r). We recall 
that, for a given eigenvalue, we recover the energy by the expressions 



F(v) 



mm s 

s>0 L 



•«/»] 



mm 

r>0 



#[/J( r ) + w /(r) 



It follows that if we scale and shift /(r) to Af(r/b) + B, a, b > 0, then the new / and if-functions are given by 

f(s) -^Af{b 2 s) + B and K^(r) — K^Mr/b). 

b z 



(26) 



The T^T- function has strong invariance since it must be employed with a changed f(r) and the impacts of the parameters 
A and B are not overlooked. For pure powers and the log potential, the if-function is completely invariant with 
respect to this family of potential-shape transformations, including scale. We shall show in the next section that if 
f(r) = g(h(r)), and g is convex, then the approximation leads to energy lower bounds; it also removes 

/ from K and suggests that inversion F — > / might be possible via Eq. (fT8|) . When g is concave, we obtain upper 
bounds. 
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III. SMOOTH TRANSFORMATIONS AND ENVELOPE APPROXIMATIONS 



In this section we consider potential shapes f(r) that may be written as smooth transformations f(r) — g{h(r)) 
of a 'basis potential' h(r). The idea is that we know the spectrum of —A + vh(r) and we try to use this to study 
the spectrum of — A + vf(r). When the transformation function g has definite convexity (g" does not change sign), 
the kinetic-potential formalism immediately allows us to determine energy bounds. This is a consequence of Jensen's 
inequality j34j . which may be expressed in our context by the following: 



g is convex (g" > 0) 
g is concave (g" < 0) => 

More specifically, we have for the kinetic potentials 

g" > =► /» > g(h(s)); 



{il>,g(h)1>)>g{y>M)), 
9" < f(s) < g(h(s)). 



(27) 



(28) 



We can summarize these results by writing /(s) w g(h(s)) and remembering that « is an inequality whenever g has 
definite convexity. The expression of these results in terms of if-functions is even simpler, for we have 

K [f] = f- 1 o / ps (g o h)- 1 o (g o h) = hT x o h = K [h] . (29) 

Thus KW « KM is the approximation we sought, that no longer depends on /. The corresponding energy bounds 
are provided by 



E = F(v) 



min { s 

s>0 L 



vg (£(*))}= min {ifW (r) + v/(r)j . 



(30) 



Many examples and results derived from Eq. (I30[) have been discussed in earlier papers [l9- 24] . It remains here for us 
to describe another equivalent but geometrical approach which originally led, without Jensen's inequality, to the term 
'envelope method'. This alternative is particularly useful when, as for the Dirac equation (35l.l36|. the kinetic-potential 
apparatus turns out to be complicated. 

We study potentials of the form /(r) = g(h(r)) and we suppose that an eigenvalue of the operator —A 
known and is given by H(v). We now consider the 'tangential potentials' given by 



where the coefficients arc 



f^{r) = a{t)h(r) + b{tl 



o(t) = g\h{t)), and b{t) = g{h{t)) - h(t)g'(h(t)). 



vh(r) is 
(31) 



(32) 

Thus the original potential f(r) is the envelope of the family of tangential potentials {f^'(r)}. Let us suppose, for 
definiteness, that g(h) is convex. Then each tangential potential is of the form ah(r) + b lying below f(r) and whose 
eigenvalues are known in terms of H(v). Thus, the eigenvalue F(v) of —A + vf(r) is bounded below by that of the 
'best' lower tangential potential, the envelope of the family of lower energy curves. Explicitly, we have 



F(v) > max[H(a(t)v) 
t>o 



b(t)i 



(33) 



If we use the expressions for <z(t) and b(t) given by Eq. (f3"2"j). then we find that the critical point satisfies H'lu) = h, 
where u — vg'(h). Meanwhile, the critical value is given by H (u) — uH'{u) -I- vg(H' (u)). But this is exactly what we 
obtain from the expressions s = H(u) — uH' '(u), h(s) = H'(u), and f(s) ps g(h(s)), that is to say 



F(v) Pa min [H{u) - uH'(u) + vg (H'(u))] . 
This corresponds to Eq. (fTTjh with the envelope approximation for the kinetic potential. 



(34) 



IV. THE ENVELOPE INVERSION SEQUENCE 



We suppose that an eigenvalue E of H — —A + vf(r) is known as function E = F(v) of the coupling parameter 
v > 0. In some cases, such as the square well, the discrete eigenvalue may exist only for sufficiently large coupling, 
v > v\. The kinetic potential /(s) may be obtained by inverting the Legendre transformation Eq. ([9]). Thus 



F(v) = min \s 
s>0 L 



vf(s)] -> f(s) = max 

V>V± 



F(v) 



(35) 
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We shall also need to invert the relation Eq. (|18l) between and by means of 

K(r) = max [F(v) — vf(r)] . 



(36) 



We begin with a seed potential p°'(r) from which we generate a sequence {f^(r)}^L of improving potential approx- 
imations. The idea behind this sequence is that we search for a transformation g so that g(f^(r)) is close to f(r) in 
the sense that the eigenvalue generated is close to F(v). The envelope approximation is used at each stage. The best 
transformation at stage n is given by using the current potential approximation /'"'(r) as an envelope basis. We 
have: 



/ = 5 l»Jo/l«J =>gW = fop 



Thus 



f[n+l] _ g[n] Q j[n] _ j Q j^\n] _ 

The resulting inversion algorithm may be summarized by the following: 
inversion algorithm 

fin] ( r ) pin] r v ) ^ R [n] ( r ) = max [pin] ( M ) _ u f\n] ^ 

U>Vi L 



f [n+l] {r) 



max 

V~>V\ 



F(v) KN( r ) 



(37) 
(38) 



The step /[™](r) -> F^(v) is effected by solving (-A + vf^)t/j = Eip numerically for E = 



V. THE INVERSION OF SPECTRAL DATA FROM PURE POWERS AND THE LOG POTENTIAL 



Suppose that the given spectral function F(v) derives from a pure-power potential shape f(r) — sgn(g)r 9 . Then 
according to Eq. (|20|) we know implicitly that the kinetic potential has the form f(s) = A(q) sgn(g) s~ q l 2 . If we effect 
the inversion by starting, for example, with a Coulomb seed /'°'(r) = —1/r, we have K^(r) = P(— l) 2 /r 2 . It follows 
from the inversion algorithm that 

/N(r) = f(K^(r)) = {p ^ ))q sgn( g ) r". 

Hence, by scaling Eq. (j2"6")l . we find ifW(r) = P(q) 2 /r 2 , and, in one more step, we have 

fW(r) = f(KW(r)) = sgn(q)r«=f(r). 



Similarly, if F(v) is generated by the log potential /(r) = ln(r), we have F(v) — vF(l) — ~v ln(u). In this case, 
f(s) = F(l) - |ln(2es). We again choose a Coulomb seed fW(r) = -1/r, for which K^(r) = P(-l) 2 /r 2 . It follows 
from the inversion algorithm that 

/W(r) = /(^°l(r)) = ^(l)-iln(2eP(-l) 2 /r 2 ) 

= F(l) - iln(2eP(-l) 2 ) +ln(r). 

By scaling Eq. flM]) we see that K^(r) = K^(r) = P(0) 2 /r 2 , and it follows that 

/ [2l (r)=/(A- [1] (r))=Mr). 
Any pure-power (or the log) seed would yield the same exact inversion in two steps for these cases. 
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VI. INVERSION OF HULTHEN SPECTRAL FUNCTIONS 



The Hulthen potential has shape f(r) = — l/(e r — 1) and s-state eigenvalues F n (v) given in Eq. Q. This potential 
has a Coulomb-like singularity at the origin as evidenced by the Coulomb- like behaviour F n (v) ~ — cv 2 of the spectral 
curves for large v. We therefore use a Coulomb seed /'°'(r) = —1/r to begin the inversion of F(v) — Fi(v). As before, 
we have from this seed that K^(r) = P(— l) 2 /r 2 . The first inversion step can be made analytically. By using f(s) 
from Eq. (|16l) with n = lwe obtain 



/W(r) = f(KM(r)) 



1 




' 4 +1 " 


1 
2 








- 1 


2 









(39) 



This step is shown in Fig. (1) as the curve clearly between the seed and the 'goal'. The iteration is repeated twice 
more numerically and the graph shows f^(r), which almost coincides with the goal f(r). We have found similar 
convergence when inverting excited-state spectral curves F n (v) for n > 1 by the same algorithm. The convergence 
is better in these cases if the seed A-function K^(r) = {P n t{—l)/r) 2 corresponds to the same radial excitation 
Pfu(— 1) = n + 1 = n as that of the energy curve that is being inverted; with this choice, the first iteration f^(r) is 
exactly given by Eq. (|39[) for all n > 1. 




FIG. 1: The inversion of Hulthen spectral data F(v) = —{v— l) 2 /4, v > 1, shown as the upper inset curve; the lower inset curve 
is the Coulomb spectral function — u 2 /4. In the main graphs the Coulomb seed = —1/r is used to start the inversion 

iteration sequence. The diagram also shows the first iteration /^'(r) and the third iteration /' 3 '(r), which almost coincides 
with the Hulthen goal ,f(r) = -l/(e r - 1). 
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FIG. 2: The ground-state eigenvalue curves F(v) for the operator H = —A + v(— 1/r + ±w(r)), where w(r) is respectively r, 
r 2 , and ln(r). The inversion of these spectral curves are shown respectively in Figs. (3 - 5). 

VII. INVERSION OF SPECTRAL DATA FROM POTENTIALS OF THE FORM f(r) = -a/r + bw(r), 

WHERE w{r) = sgn(q)r q OR ln(r). 

We now consider three examples of spectral data arising from singular potentials of the form 

f(r) = \-bw(r), wherew(r) = r, r 2 , ln(r). (40) 

r 

Graphs of the corresponding ground-state energy functions F(v) for a — 1 and b = | are shown in Fig. @. The 
inversion algorithm of section (4) has been applied to each of these spectral curves and the results are shown in 
Figs. (j3H5j) - In each case the seed potential was the Coulomb potentinal /'°'(r) = —1/r. The curve lying above the 
goal is the first iteration f^(r); the third iteration is also shown and is in each case very close to the goal of the 
example, namely f(r) itself. 

VIII. UNIQUENESS 

An important issue concerning the inversion F — > f is whether there is more than one inverse / for a given F. 
We tackled this question earlier, in Ref. [lij §6, for bounded symmetric potentials in d = 1 dimension that were 
monotone non-decreasing on the half-line x > 0: under these conditions, the inverse was shown to be unique. It is 
straightforward to prove a similar theorem for bounded potentials in three dimensions. However, a new theorem is 
needed for singular potentials. We suppose in d = 3 dimensions that the potential shape f(r) is attractive, spherically 
symmetric, and Coulombic for small r; more specifically, weshall assume that /(r) can be written in the form 

/W = — , where g(0) < 0, g'(r) > 0, (41) 
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5 
4 




r 



FIG. 3: The inversion of the spectral data shown in Fig. (2) for the Coulomb-plus-linear potential f(r) = —l/r+ hr. The 
diagram shows the seed /'°'(V) = — 1/r, the first iteration f^(r), and the third iteration /' 3 '(r), which almost coincides with 
the goal, the exact potential f(r). 



and g(r) is not constant. Examples of this class of singular potential shapes are Yukawa g(r) = — e~ ar , Hulthen 
g(r) — —r/(e ar — 1), and linear-plus- Coulomb g(r) — —a + br 2 , with a, b > 0. We consider the ground state ip(r,v) 
which we suppose to be a normalized radial function so that 

y* oo 

/ i/j 2 (r,v)r 2 dr = 1. (42) 
Jo 

With these assumptions, we shall prove the following 

Theorem 1 The potential shape f(r) in H = —A + vf(r) is uniquely determined by the ground-state energy function 
E = F(v). 

Proof 

For this singular class of potentials it is helpful to consider an additional probability density on [0, oo) of the form 
cip 2 (r 7 v)r and to define the integrals 

/>a />oo 

q(a, v) = / ijj 2 (r,v)rdr and I(v) = / ip 2 (r,v)rdr. (43) 
Jo Jo 

In this sense, Q(a,v) = q(a,v)/I(v) is the probability mass on [0,a\. Our hrst task is to establish a concentration 
lemma to the effect that, for each fixed a, this probability Q(a, v) — > I as v increases without limit. From Eq. (0) we 
have 

/>oo y* oo 

F'(v)= *p 2 (r,v)r 2 f(r)dr = ^ 2 {r )V )rg{r)dr > l(v)g(0) (44) 
Jo Jo 
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2 



f(r) 



-3 



-4 




FIG. 4: The inversion of the spectral data shown in Fig. (2) for the Coulomb-plus-oscillator potential f(r) = — 1/r + \ r ' 2 ■ The 
diagram shows the seed = the first iteration f^(r), and the third iteration /' 3 '(r), which almost coincides with 

the goal, the exact potential f(r). 



and 



Thus we find for each a such that g{a) > g(0), 



F'(v) = / tp 2 (r,v)rg(r)dr + / ip 2 (r,v)rg(r)dr > q(a,v)g(0) + (I(v) — q(a,v))g(a). 



,^n< , •> 9(a) - F'(v)/I(v) 

1 > Q(a,v) = — . > 



I(v) 



g(a) - g(0) 



(45) 



(46) 



For large v the potential shape is spectrally dominated by the Coulomb term vf(r) ~ vg(0)/r. The corresponding 
wave function, eigenvalue, and I(v) are given asymptotically by 



i/}{r,v) ~ce V9{0)r/2 , F(v) 



and I( V )J^M. 



(47) 



It follows that F'(v)/I(v) ~ g(0). By using these asymptotic expressions in Eq. (|4U)) and Eq. (l4"4")l we arrive at the 
limit 



lim . 

t)— s-oc \ I(v) 



3(0). 



Thus, for each fixed a such that g(a) > 5(0), we have the radial concentration lemma 

,^ n , v ^ g(°) - F'(y)/I(v) 

1 > Q(a,u) > -r >■ 1. 

9W - .9(0) 



(48) 



(49) 
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f(r) 



-2 



-4 




FIG. 5: The inversion of the spectral data shown in Fig. (2) for the Coulomb-plus-log potential f(r) — —1/r+i ln(r). The 
diagram shows the seed = — 1/r, the first iteration /'^(r), and the third iteration f^(r), which almost coincides with 

the goal, the exact potential f(r). 



Now let us suppose that the two potential shapes g{r)/r and (g(r) + 7(r))/r both have the same ground-state energy 

function F(v). We shall prove that 7 = 0. By considering large v we see from Eq. (|47|) that 7(0) = 0. We consider 
the two Hamiltonians H = —A + vg(r)/r and Hi = —A + v(g(r) + 7(r))/r with respective lowest normalized radial 
eigenstates if){r, v) and tf>(r, v), J(v) = J °° 4> 2 (r, v)rdr, and the common eigenvalue curve is F(v). From the equalities 



and the variational inequalities 



we deduce the complimentary inequalities 



(^,H^)>F(v) < (<i>,H<P) 



(4>,"f(r)/r<f>) = / (j) (r,v)^/(r)rdr < < / ip (r,v)j(r)rdr = (ip,"f(r)/r ip). 
Jo Jo 



Thus 



— — — / <ft (r,v)^(r)rdr < < —. \ ip (r,v)~f(r)rdr 
J{v)Jo I \ v ) Jo 



(50) 



(51) 



(52) 



(53) 



By a sweeping process starting at r = we now show that Eq. (|53|) implies that 7 = 0. We know that 7(0) = 0. 
Clearly, if 7(7*) does not change sign, then it must vanish identically. Let us suppose that 7(7") first changes sign at 
r = a; more specifically (and without loss of generality), we suppose 7(f) > for r £ (0, a), and that 7(a) = 0. We 
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first note that if g(a) — g(0), g(r) would be constant on [0, a] and 7(1") would have to vanish on [0, a] in order for 
idi 7 ") + l{ T ))/ r to remain in the same potential class. Thus we understand r — a to be the first point where 7(7") is 
zero, after it has been positive. We now increase v so that both of the following inequalities are satisfied 

F '( v ) , x F '( v ) f \ / \ / \ 

In view of the radial concentration lemma, further increase in v will force both integrals in Eq. (|53p to be positive. 
This cannot be unless 7 = 0. This establishes Theorem 1. D 



IX. CONCLUSION 



We are presented with a spectral curve, F(v), which shows how a discrete eigenvalue of the Hamiltonian H = 
— A + vf(r) depends on the coupling parameter v > 0. From this data, we reconstruct the underlying potential shape 
f(r). This 'geometric spectral inversion' is effected by a functional sequence in which the unknown potential shape 
/(r) is first regarded as a smooth transformation g (/^(f)) °f a seed potential /'°'(r). The paper focuses on problems 
that involve singular potentials for which the Coulomb seed /'°'(r) = — 1/r is very effective. The method of potential 
envelopes is used to generate a sequence of approximations {f^(r)}^L for which the corresponding transformations 
gradually approach the identity. At this stage of the investigation we offer an explicit heuristic inversion algorithm 
rather than an abstract convergence theorem. However, for the class of attractive central potentials with shapes 
f(r) = g(r)/r, with g(0) < and g'(r) > 0, we prove that the ground-state energy curve F(v) determines f(r) 
uniquely. 
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